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Abstract 

We present a new approximation technique for quantum field the- 
ory. The standard one-loop result is used as a seed for a recursive 
formula that gives a sequence of improved Gaussian approximations 
for the generating functional. In a different setting, the basic idea 
of this recursive scheme is used in the second part of the paper to 
substantialy speed up the standard Monte Carlo algorithm. 
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1 Introduction 



Quantum field theory is compactly written in terms of path integrals. Path 
integrals have been specially usefull in dealing with symmetries, quantizing 
gauge theories, etc. At some point, however, we need to calculate the path 
integrals. The problem is that the only path integrals that we know how to 
solve correspond to free theories (Gaussian integrals). We would be in a very 
sorry state if we didn't have a generic approximation scheme at our disposal. 
Semi-classical or loop expansion is just such an approximation scheme. In 
fact, much of what we know about quantum field theory comes from one- 
loop results. The one-loop result is obtained by Taylor expanding the action 
around classical fields and disregarding cubic and higher terms. In this way 
the path integral is approximated by a Gaussian. 

In this paper we will consider another Gaussian approximation to the 
path integral. Unlike the one-loop result, here we will Taylor expand about 
the average field (p = (0). We shall show that this leads to an improved 
approximation given in terms of a recursive relation. 

2 The Gaussian Approximation 

The central object in quantum field theory is the generating functional Z[J]. 
Functional derivatives of Z[J] with respect to the external fields J(x) give 
the Green's functions of the theory. The generating functional is determined 
from the (Euclidian) action S[(/>] through the path integral 



where d is the dimension of space-time. We are interested in looking at a set 
of approximations to the above path integral. The approximations are valid 
in all d. In this letter, however, we are going to look at the simpler case of 
d — theories, where it is easy to compare our results with exact numerical 
calculations. In d = functionals become functions, and the path integral 
reverts to a single definite integral over the whole real line. 




(1) 



The integration measure is, formaly, simply 




(2) 




(3) 
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An even more useful object is W(J) — the generator of connected diagrams, 
defined by 

Z(J) = Z(0) e -* w{J) . (4) 

In statistical mechanics parlance this is the free energy. The quantum average 
of the field is 

*=(<!>) = — HjW{J). (5) 

In the Gaussian approximation, we Taylor expand the action in the path 
integral around some reference point re f, and keep terms that are at most 
quadratic in — ref . Thus, we use 

S(0) « S(0 ref ) + S'(0 ref ) (0 - ref ) + ~ S"(0 ref ) (0 - rcf ) 2 . (6) 

The integral in (JHJ) is now a Gaussian and we find, up to an unimportant 
constant, that 

WW J, ref ) = S(0 rcf ) - J0 rcf + 5 In S"'(0 ref ) - \ {S ' { t c 'l " • (7) 

2 2 6"(0 re f) 

For this approximation to make sense, the integral must get its dominant con- 
tribution from the vicinity of the reference point re f. The standard Gaussian 
approximation corresponds to the choice re f = c i aS s (J), where c i aS s is the 
solution of the classical equation of motion S' = J. The classical solution is 
the maximum of the integrand in (JHJ). This specific choice of re f gives us the 
standard one-loop result 

W(J) « Wi(J) = 5(0 class ) - J0 class + ^ ln5"(0 class ) . (8) 

As is well known, loop expansion is just an expansion in powers of h. The 
one-loop result gives us the first quantum correction to classical physics. 
From now on we set h = 1 . 



3 Improving the Gaussian Approximation 

In this section we will choose a different expansion point re f for our general 
Gaussian formula (j2J). The idea is to expand around the average field (p. Al- 
though the classical solution gives the maximum of the integrand, expansion 
around (p gives a better approximation for the area under the curve. This is 
in particular true for large values of J. We will work with 4 theory in d = 0, 
whose action is given by 

S(0) = ^0 2 + ^0 4 • (9) 
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The classical equation of motion is now a cubic algebraic equation. We easily 
find the unique real solution. In this way we get a closed form expression 
for the one-loop approximation W\{J). The Gaussian approximation around 
the average field ip is simply 

S(v) - Jv + \^ (v) -\^-f. ( 10 ) 

To be able to calculate this in closed form we need to know <p{J), which 
is tantamount to knowing how to do the theory exactly, since (p and its 
derivatives give all the connected Green's functions. The use of equation 
(JTUJ) comes about when one solves it iteratively. We use equations © and 
(J7|) as the basis for the following iterative process 

<Pn+l(J) =~Jj U <:-»(•/• <Pn(J)) ■ (11) 

Differentiating (jJJ) we find 

gQgn) - J _ 

S"( Vn ) 

S"'(p n ) ( fS>(tp n ) -J\\ _J_] *P* . (12) 



V^n+l — <-fr 



S"(<p n ) J S»{<p n ) J dJ ■ 

For the seed of this iteration we choose the classical field, i.e. <po = 0ciass- 
In this way we obtain a sequence of points ipo, ipi, <p2, ■ ■ ■ or equivallently of 
approximations to the connected generating functional Wi, W 2 , W 3 , . . . given 
by W n+ i(J) = Wq^s^J, (p n (J)). The idea behind this is obvious — we want 
to obtain a sequence of W n (J)'s that give better and better approximations 
to W(J). The following two figures show that this realy works. 

From Figure [T] we see that the sequence of <£> n 's converges to <poo 0. 
Note that ip^ ^ (p. The reason for this is obvious: We used the Gaus- 
sian approximation WGauss i n defining our recursive relation, and there is 
no reason to expect that this converges to the exact result. It does, how- 
ever, converge and (p^ represents an excellent approximation to p>. To see 
this, in Figure El we have ploted the ratio | \ . This ratio represents 
a direct measure of the improvement of approximations in going from the 
one-loop result W\(J) = Wq^^J, <po(J)) to our improved Gaussian result 

w 00 (j) = w G ^M^oo(J))- 

In iterating (|12fl we necessarily discretize the J's. The coursness of this discretization 
effects the speed of convergence of the (p n 's. The standard way out of this problem is to 
introduce a small mixing papameter e. Instead of the recursive relation y n +i = f{tp n ) one 
then consideres ip n +i = (1 — e)<-Pn + e f( l fn)- 
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Figure 2: Plot of the ratio \((p — <p ) / (<p — <^oo)| as a function of J. This is 
a direct measure of how the new approximation outperforms the standard 
one- loop result. Here we had g — 1. 



As we have already mentioned, our new approximation was tailored to 
work well for large J. This can be read off directly from Figure |2] For 
example, for J = 15 the improved Gaussian is about one hundred times 
better than the one- loop result. The new approximation is poorest for J w 2, 
but even there it beats the old approximation by a factor of seven. Most of 
the time we are interested in working with small or zero external fields. In 
the vicinity of J = the new approximation is fourteen times better than 
the old one. 
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4 Monte Carlo 



The aim of our investigations so far has been to develop better analytic ap- 
proximation schemes that can be applied to general quantum field theories. 
We worked in d = in order to be able to make a simple comparison with 
exact (numerical) results. In this section we will look at the numerical tech- 
niques themselves. We used the Monte Carlo algorithm pQ for calculating 
path integrals. In d = Monte Carlo is not the most efficient way to do 
things — its advantages become apparent as we look at larger and larger 
numbers of integrations. We use Monte Carlo in order to investigate the 
algorithm itself in light of what we have learned in the previous two sections. 

We start with a brief introduction of the method. In order to calculate the 
definite integral J f((p)d<p we choose a non-negative function p(4>) normalized 
so that J p(<f>) d(j) — 1. Therefore, p((fi) is a probability distribution. The 
integral is now 

md(p= I m pWd ^(i) P ' (13) 

where (F) p represents the mean value of F with respect to the probability 
distribution p. Therefore, the integral of / is given as the mean value of f /p 
on a sample of random numbers whose probability distribution is given by 
p. In practice, this mean value is estimated using a finite number N mc of 
Monte Carlo samples, and the error of such an estimate is itself estimated to 
be <jf/p = \J<2 2 f/ p i where the variance equals 

((£\ 2 \ .It 2 

a 9 



f/P N - I 

1 v mc x 



l v 

(14) 



The central limit theorem guarantees that the Monte Carlo algorithm con- 
verges to J f d(p for an arbitrary choice of distribution p. The only condition 
that must be met is a^, p < oo. This freedom in the choice of p is used 
to speed-up the convergence of the algorithm. The speed of convergence is 
measured by the efficiency £ , given by 



1 a f/ P 



where T represents the total computation time. Note that a hundred fold 
increase of efficiency corresponds to one extra significant figure in the final 
result. 
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In our calculation we chose p(4>) to be the Gaussian normal distribution 



where a and a completely determine the distribution. There are two reasons 
for using this distribution. First of all, the function we are integrating can be 
approximated by a Gaussian over a wide range of parameters J and g. A good 
choice of a and a makes f/p almost constant over the range of integration, 
thus making the variance small. Second, there exists a specific algorithm 
for generating random numbers conforming to a Gaussian distribution. The 
Box-Muller algorithm [2j is much more efficient than the standard Metropo- 
lis algorithm [Sj since it doesn't give rise to autocorrelations of generated 
numbers. In the Metropolis algorithm autocorrelations can be pronounced, 
and their removal substantialy slows down the simulation. 

The choice of probability distribution has a great effect on the efficiency. 
For example, the efficiency corresponding to the uniform distribution on the 
interval <p £ [ — 100, 100] is 3.5 10 10 times smaller than the efficiency achieved 
by the Gaussian distribution centered at a = c i ass with optimal choice of 
width a. Having chosen p to be a Gaussian, the computation time T depends 
only on the number of Monte Carlo samples iV mc . Therefore, in our case, 
maximalization of efficiency is equivallent to a minimalization of the variance 

a )lv 

In the previous sections we saw that it is even better to expand around 
if. In the Monte Carlo setting this should translate into a further increase 
in efficiency. This is precisely what we see. By varying the center of the 
Gaussian a (always using optimal width for that given a), we find maximum 
efficiency precisely at a = (f as we can see in Figure 01 

Figure |U compares the efficiencies Eq of simulations about c i aS s and £q 
for simulations about if for various values of J. It is seen that we get a 
two fold improvement in efficiency. This may not seem spectacular, and in 
d = it realy is not. However, once we consider theories in d > we are 
dealing with true path integrals. If we approximate the path integral with N 
integrals then the expansion around f gives a jump in efficiency of 2^. Even 
for a modest simulation with N = 20 this corresponds to an increase of six 
orders of magnitude. 

The problem with this calculation is that we already need to have the 
exact result for <f in order to get the stated increase in speed. The way out 
is obvious and is reminiscent of the step we made in the previous section in 
going from (fTU|) to (jllj) . Therefore, we need to start Monte Carlo with a 
Gaussian distribution centered about c i aS s- After a while this gives us an 




(16) 
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Figure 3: The variance as a function a. The plot is for g — 10, J = 1. The 
variance is minimized for <p(l) = 0.376799 (black dot). The classical field is 
0ciass(l) = 0.614072 (grey dot). 
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Figure 4: The ratio £q/£c as a function of J for g = 1. 



approximation to cp, say (p™ c . Using this as the center of a new probabil- 
ity distribution we obtain (/9™ c , etc. Unlike the series <po, (fx, y?2, • • • of the 
previous section, this one necessarily converges to the exact result — even 
ordinary Monte Carlo does that. The improved Monte Carlo scheme, how- 
ever, can be tailored to yield an efficiency very near to the ideal value £q 

i 
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5 Conclusion 



We have looked at two different ways how one can take advantage of the 
(rather intuitive) fact that in quantum field theory Gaussians are best cen- 
tered about the average field ip = ((f)). We cast the Gaussian approximation 
about ip as a recursive relation. Working on 4 theory in d = we have shown 
that the iterates of this equation present better and better approximations 
to W(J). This sequence of approximations ends with Woo(J), i.e. the best 
Gaussian approximation. The first iterate in this sequence is the standard 
one-loop result. The second (W2) is not much more complicated, and is al- 
ready much better than the one-loop approximation. In looking at theories 
in d > it may be possible to use W 2 to get a better analytic approximation 
to (say) the effective potential. 

The second use for the newly centered Gaussians is in Monte Carlo cal- 
culations. The Monte Carlo algorithm is most efficient when one generates 
random numbers through Gaussian probability distributions centered about 
if. We have shown that this can substantialy speed up the algorithm. For 
an iV-fold integral the speed up is roughly 2 N . We are currently working on 
applying the improved Monte Carlo algorithm to models in d > 1. 
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